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The relative stability of different high-pressure phases of various Cesium Halides is studied from 
first principles and analyzed using the Landau theory of phase transitions. We present results for 
Csl, CsBr, and CsCl up to pressures of 100 GPa. A cubic-to-orthorhombic transition, driven 
by the softening of an acoustic phonon at the M point of the Brillouin zone, is competing with 
the cubic-to-tetragonal martensitic transition typical of these compounds. The phonon softening 
takes place only in Csl and CsBr at a residual volume of V/ Vb=0.64, 0.52 respectively. A cubic- 
to-tetragonal instability is found instead to occur at V/Vo » 0.54 for all the compounds considered 
here. The orthorhombic phase is stable only in Csl, whereas it is taken over by the tetragonal one 
in the case of CsBr. Our analysis reveals the essential role played by the phonon-strain coupling in 
stabilizing the orthorhombic phase and in making the corresponding transition first-order. 



I. INTRODUCTION 

Since the pioneering work of Madelung and Ewald, |l]] 
a huge amount of experimental and theoretical work has 
been devoted to alkali halides which are the simplest and 
most representative ionic solids. The interest in the high- 
pressure properties of these materials has been recently 
renewed by the introduction of the diamond anvil cell 
technology [|| which has opened the way to exploring pre- 
viously unaccessible portions of their phase diagram. Csl 
is the softest among alkali halides and it has the smallest 
optical gap. For these reasons it is an ideal candidate 
to display a pressure-induced band-overlap metallization 
which is actually found at a pressure P rj 110 GPa. 
H Looking for this effect, new and unexpected crystal 
phases have been found to be stable in the pressure range 
of a few tens GPa. 

The first systematic studies of the high-pressure phases 
of Cesium halides date back to 1984. From energy- 
dispersive X-ray diffraction experiments Huang et al. Q 
and Knittle et al. || were able to observe that Csl and 
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CsBr (whose low-pressure stable phase is the cubic B2), 
undergo a transition which lowers their cubic symme- 
try when an applied pressure reduces their residual vol- 
ume to v — V/Vq ps 0.54, Vq being the equilibrium vol- 
ume. The observed transition pressures are of ps 39 and 
53 GPa for Csl and CsBr respectively. || A similar insta- 
bility has been also reported in CsCl, at about 65 GPa. 
H This symmetry lowering manifests itself through the 
splitting of the (110) X-ray diffraction line into two peaks 
which arc interpreted as corresponding to the (101) and 
(110) incquivalent reciprocal-lattice vectors of a tetrag- 
onal structure. These authors did not observe any de- 
tectable discontinuity in any physical observable, and 
they conjectured therefore that the phase transition was 
of second order. However, general group-theoretical con- 
siderations indicate that a transition from a cubic to a 
tetragonal phase in which the only order parameter is the 
unit cell shape, i.e. strain, cannot be of second order. || 
We conclude therefore that a weak first-order character 
of the transition must have escaped the analysis of the 
experimental data. Claims that this tetragonal phase 
would undergo a further transition to an orthorhombic 
structure Q| (space group D\ h ) were not confirmed by 
subsequent theoretical work (|-[l(J . 

The first-order character of the cubic-to-tetragonal 
transition in Cesium halides was confirmed by the very 
first theoretical investigations performed by Vohra et al. 
Q within a Born-Meyer semiempirical approach. Ac- 
cording to these authors, the instability of the CsCl crys- 
tal structure is related to the softening of the shear con- 
stant, c s = \(cn — C12), due to the competition between 
the long-range Coulomb attraction and the short-range 
repulsive energy which becomes more and more impor- 
tant as the pressure increases. They also found that the 
tetragonal structure is stable against further orthorhom- 
bic distortions of symmetry corresponding to the space 
group D\ h \ i.e., for any given value of a and c/a, the min- 
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imum of the lattice energy was always found abb/ a = 1 in 
the whole pressure range explored. Some investigations 
on Csl have also been carried out from first principles. 
The work of Christensen and Satpathy || — performed 
within the Linear Muffin-Tin Orbital approximation — 
and that of Baroni and Giannozzi p0| — in which ab-initio 
pseudopotentials were used — confirmed the overall pic- 
ture proposed by Vohra et al. 

In 1989, new and unexpected features emerged from 
the X-ray diffraction experiments performed by Mao et 
al. |ll[] on Csl. Using a sophisticated experimental set-up 
based on a synchrotron radiation source, they were able 
to identify a distortion of Csl from the cubic B2 to the 
hexagonal closed-packed structure passing through an in- 
termediate orthorhombic phase whose assigned symme- 
try (C^,), however is different from that guessed in Ref. 

0. It was observed that the (110) peak splits into three 
components, whereas two additional lines appear on the 
low- and high-energy sides of the triplet, corresponding 
to two diffraction peaks which are forbidden in the B2 
structure. This diffraction pattern can be explained by 
assuming a 2 x 1 reconstruction of the cubic cell character- 
ized by the gliding of one (110) plane, as depicted in Fig. 

1. The onset of the phase transition, characterized by the 
broadening of the diffraction lines of the cubic structure, 
is observed at pressures as low as 15 GPa. The observed 
splittings and peak intensities evolve continuously from 
15 to 100 GPa, showing a single orthorhombic phase with 
variable parameters. However, below 45 GPa the split- 
tings are not large enough to allow an unequivocal iden- 
tification of the orthorhombic structure The transition 
from the cubic to the orthorhombic phase was observed 
to be continuous within the experimental resolution. 

In a recent paper | p2[ we have shown that the distortion 
leading to the orthorhombic structure observed by Mao 
et al. [ pT| is driven by the softening of an phonon at 
the M point of the Brillouin zone (BZ). This softening is 
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related to the ferroelastic instability responsible for the 
previously assumed cubic-to-tetragonal transition. The 
shear constant is in fact proportional to the square of 
the sound velocity along the (110) direction for vibrations 
polarized along (110). This observation suggests that a 
large reduction of the magnitude of the shear constant 
could favour the softening of a transverse phonon along 
the (110) direction. In fact, the gliding of the (110) plane 
depicted in Fig. 1 represents a lattice distortion which is 
very similar (but not strictly equal) to one of the two de- 
generate ionic displacement patterns of the M§ acoustic 
phonon mode. We are thus led to identify the amplitude 
of this mode as the order parameter of the transition, in 
Landau's sense. 

In this paper we present a complete account of our pre- 
vious work on the high-pressure phases of Csl, [JL2| and 
its extension to other Cesium halides, CsBr and CsCl. 
Our approach is based on a combined use of the Landau 
theory of phase transitions to classify the possible low- 
symmetry phases and to determine the form of the inter- 
actions which can lead to them, and on first-principles 
calculations of the relevant interaction constants and en- 
ergy differences as functions of the applied pressure. Our 
calculations are based on density-functional theory to de- 
termine crystal energies and their first derivatives with 
respect to atomic displacements (forces) and macroscopic 
strain (stress), and on density- functional perturbation 
theory to determine the second derivatives (essentially, 
the phonon frequencies). In Sec. II we present our the- 
oretical framework and computational details. Sec. Ill 
contains our results. Sec. IV contains the conclusions. 

II. THEORY 
A. Landau theory of phase transitions 

In the Landau theory of phase transitions, the rele- 
vant thermodynamical potential, J 7 , describing the rela- 
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tive stability of two phases is expressed as a power series 
of the so called order parameter, u, whose value is dif- 
ferent from zero in the low-symmetry (high-pressure, in 
the present case) phase, while u vanishes in the high- 
symmetry (low-pressure) phase: 

T(u) = JF + A 2 u 2 + A 3 u 3 + A 4 m 4 + • • • , (1) 

where the linear term in u is missing to ensure the sta- 
bility of the low-pressure phase. Let us suppose for the 
moment that the fourth-order coefficient, A4, is posi- 
tive so that it can ensure the global stability of the 
system, and let us neglect all higher-order terms. If 
A 2 > 0, the relative stability of the low- and high- 
symmetry phases is determined by the sign of the dis- 
criminant A 6 = A3 2 - AA 2 A i . If A 6 < the high- 
symmetry (u = 0) phase is stable, while for A 6 > 
Eq. (1) has a minimum for u ^ and the low-symmetry 
phase prevails. In this case, the transition from the high- 
to the low-symmetry phase is discontinuous, and it is said 
to be first-order. Continuous second-order transitions are 
possible if, because of symmetry, the third-order coeffi- 
cient vanishes identically: A 3 = 0. In this case the high- 
or low-symmetry phases are stable according to whether 
the second-order coefficient, A 2 , is positive or negative. 
Suppose now that A4 is negative in the region of the 
phase diagram where A 2 softens and that A 5 = 0. In this 
case, the global stability of the system must be ensured 
by the sixth- and higher-order terms and a first-order 
transition to the low-symmetry phase would then occur 
when A 8 = A 4 2 - AA 2 A 6 = 0. 

For structural transitions where the order parameter 
is the amplitude of a lattice distortion, A 2 is essentially 
the square of the phonon frequency associated to that 
distortion. In the present case, the phonon which goes 
soft is the Mg" mode at the M point of the BZ. The small 
group of the M point — whose coordinates are (|^0) — is 
Z?4/j and its star is made of three equivalent points. 
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transforms according to a doubly degenerate irreducible 
representation of D^) ll so that the order parameter asso- 
ciated with this mode is six-dimensional: u = {uj}i=i,6' 
The M 5 ~ representation is odd with respect to inversion, 
so that no third-order (or, more generally, odd-order) 
terms must be present in Eq. (1) in order to ensure in- 
variance with respect to the (cubic) symmetry group of 
the low-pressure phase. Following Landau's argument, 
we conclude that A3 = 0, and A$ = 0. Assuming that 
the sixth- and higher-order terms in Eq. (1) are positive 
in the neighborhood of the transition, this is then first- 
or second-order, according to whether the fourth-order 
coefficient is negative or positive, respectively. 

B. First-principles techniques 

In this paper we will ignore any entropy effects and 
restrict ourselves to zero temperature. In this condi- 
tion, the implementation of the above phenomenolog- 
ical scheme requires the knowledge of the system en- 
ergy as a function of the atomic displacements (phonon 
amplitudes) and unit-cell distortions (strain). Density- 
Functional Theory (DFT) is a very reliable and compu- 
tationally viable tool for the study of the energetics of 
simple materials. We refer the reader to the many avail- 
able review papers on this subject [[l4| and we summarize 
here the main points. The local-density approximation 
(LDA) to DFT provides a practical, though approximate, 
way to map the ground-state properties of a system of in- 
teracting electrons onto those of a system of noninteract- 
ing electrons, subject to an effective external field which 
depends self-consistently on the ground-state electronic 
density. DFT-LDA allows to calculate in a rather ac- 
curate and efficient way the energies of materials within 
the Born-Oppenheimer approximation, along with their 
first derivatives with respect to external parameters, such 
as e.g. nuclear coordinates (i.e. atomic forces) jL5| and 
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unit-cell shape and volume (i.e. stress). JTsJ By lineariz- 
ing the DFT-LDA equations with respect to the strenght 
of an external perturbation, the second derivatives of the 
energy can be conveniently obtained by solving a suit- 
able set of linear equations. 0,|l8| This is the essence 
of the Density-Functional Perturbation Theory (DFPT) 
approach to lattice dynamics, which we will employ to 
determine the phonon spectra of the system considered, 
as functions of the crystal volume, i.e. of the applied 
pressure. 

Our calculations are performed in the framework of 
the plane-wave pseudopotential method. For Cs and 
I we have used the same pseudopotentials as in Ref. [[To) . 
Pseudopotentials for Br and CI have been generated us- 
ing the method originally proposed by von Barth and 
Car. |l9| For Cs, electrons up to Ad have been treated 
as frozen in the core. The inclusion of the Cs 5s and 5p 
states into the valence shell is necessary to obtain sensible 
results. In fact we found that a Cs pseudopotential which 
sustains only the 6s electron as a valence electron yields 
a monotonic decrease of the crystal energy at decreas- 
ing volume, with no stability at all. An explicit account 



for the non-linear core-corrections, |20| suggested by the 
large size of the Cs core, does not substantially improve 
the results. The inclusion of Cs 5s and bp orbitals in 
the valence does not increase dramatically the numerical 
burden. In fact, the spatial extension of these orbitals is 
comparable to those of the anion s and p valence wave- 
functions, so that the plane-wave basis set necessary to 
describe the latter is also adequate for the former. 

Plane waves up to a kinetic energy cutoff of 25 Ry were 
included in the basis set. The electron-gas exchange- 
correlation energy and potential used to implement LDA 
are those determined by Ceperley and Alder as in- 
terpolated by Perdew and Zunger. |2^] Brillouin-zone 
integrations have been performed using sets of special 
points corresponding to the (666) Monkhorst-Pack mesh. 
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p3| This mesh — which corresponds to 10, 12, and 24 
points in the cubic, tetragonal, and orthorombic struc- 
tures respectively — has been explicitly checked for con- 
vergence for all the relevant structural properties of in- 
terest here. 

III. RESULTS 

A. Zero-pressure properties 

The low-pressure structural properties of the materials 
investigated have been determined by fitting the calcu- 
lated crystal energies to the Murnaghan equation of state. 
Ipij DFPT allows a straightforward determination of the 
electronic contribution to the static dielectric constants 
and Born effective charges. jlTj The theoretical predic- 
tions for these data are compared with experiments in 
Table |. The resulting accuracy on the predicted equi- 
librium lattice constants and bulk moduli is of the order 
of 3% and 10% respectively, i.e. in the typical range of 
DFT-LDA calculations. Dielectric constants are system- 
atically overestimated by 20%, a well known drawback 
of the LDA. |2f| Using the theoretically determined equi- 
librium lattice parameters, the phonon dispersions of the 
three salts have been calculated by DFPT along the main 
symmetry directions. In Fig. we display the calcu- 
lated dispersion of Csl and compare them with neutron- 
diffraction data. Selected frequencies for Csl, CsBr, and 
CsCl are reported and compared with experiments in Ta- 
ble [n]. The agreement is quite satisfactory and gives us 
confidence in the predictive power of our calculations. 

B. Phonons at high pressure 

In order to substantiate the hypothesis of a mode soft- 
ening occurring upon increasing pressure along the (110) 
acoustic branch, we have calculated the corresponding 
phonon dispersion for several values of the crystal resid- 
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ual volume, v. In Fig. |3| we display our results for Csl. 
The Mg" acoustic phonon softens at v* = 0.638, corre- 
sponding to a pressure P* « 23 GPa, as a consequence 
of the incipient softening of the sound velocity of one of 
the transverse branches. The same behavior is observed 
in CsBr at v* = 0.520 while CsCl does not display any 
tendency to softening upon increasing pressure. The de- 
pendence of the Mg~ acoustic phonon frequency upon the 
residual volume is displayed in Fig. [l| for Csl, CsBr and 
CsCl. 

C. High-pressure phases 

1. Tetragonal phase 

As a first step we have reinvestigated the cubic-to- 
tetragonal transition described in the introduction. In 
Fig. ^ we display the zero-temperature enthalpies vs. 
the c/a cell parameter calculated at constant pressure 
for Csl, CsBr and CsCl. In Csl and in CsBr the exis- 
tence of a first-order transition is clearly visible with a 
second minimum occurring at c/a > 1 above some criti- 
cal pressure, Ptet , and separated from the high-simmetry 
phase by a well defined enthalpy barrier. The transi- 
tion pressures are Ptet = 44 and 58 GPa for Csl and 
CsBr respectively, with a very small volume discontinu- 
ity 0.5% and « 0.7% respectively). In Fig. || we 
report the c/a parameter vs. P/ Ptet for Csl and CsBr. 
In CsCl the stability of the low-simmetry phase cannot be 
unambiguously established within our present accuracy. 
An instability towards the tetragonal phase is signalled 
by the flatness of the enthalpy as a function of c/a for 
pressures 70 GPa. At higher pressures, however, this 
instability seems to weaken, and calculations for pres- 
sures as high as 100 GPa do not give any evidence of a 
stable tetragonal phase. 
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2. Orthorhombic phase 



As we have seen in Sec. IIA, the vanishing of the 
M§ phonon frequency signals the onset of a phase tran- 
sition whose order parameter is the phonon amplitude 
and which can be second-order because no third-order 
terms are present in the expansion of the crystal energy 
in powers of the order parameter. In order to investi- 
gate the character of the transition and determine the 
possible low-symmetry phases, it is necessary to find out 
the expression of the crystal energy in powers of u to 



the lowest meaningful order, i.e. n = 4. In Table III 
we display the wavevectors and the polarizations corre- 
sponding to one possible choice for the components of 
the six-dimensional order parameter. The corresponding 
fourth-order invariant polynomials can be obtained by 
using standard group theoretical techniques. The it's re- 
alize an irreducible representation of the symmetry group 
of the high-symmetry phase (space group 0\). The num- 
ber of fourth-order invariants is equal to the number of 
times the identical representation of Oh is contained in 
the symmetric part of [Afg~] 4 . A simple exercise shows 
that this number is equal to four. A straightforward way 
to obtain them is to symmetrize all the possible trans- 
lationally invariant fourth-order monomials with respect 
to Oh- In order to satisfy translational invariance, the 
starting monomials, UiUjU^ui must satisfy the relation: 
<li + qj + <lk + q; = G, where is the wavevector associ- 
ated to the i-th component of the order parameter, and 
G is a reciprocal-lattice vector of the cubic phase. The 



four independent invariants are summarized in Table IV 
The resulting expansion of the crystal energy reads: 



^(u) = To + A 2 ]T + A \ p 4(u) + 0(u 6 ). (2) 

t=l,6 i=l,4 



The possible low-symmetry stable phases compatible 
with this expression for the crystal energy have been 
worked out in Ref. (B5[. Following the notation of Ref. 
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|B5|, these stable phases can be classified as in Table 0, 
where each configuration corresponds to a particular di- 
rection in the six-dimensional order parameter space, i.e. 
to a particular combination of phonon displacements. In 
the specific case of Csl, the four fourth-order coefficients 
appearing in Eq. (2) have been fitted to the results of 
DFT-LDA calculations performed for lattice distortions 
along some of the order-parameter directions listed in Ta- 
ble [v|. In fact, the coefficients of the crystal energy expan- 
sion are easily obtained through total energy and stress 
calculations. With the numerical values so obtained, we 
have verified that at the crystal volume where the M§ 
phonon frequency softens, the only minimum of the en- 
ergy given by Eq. (2) occurs along the Pll direction, 
whereas directional minima along all the other directions 
are actually saddle points. This observation seems to 
suggest that a second-order transition to a phase of tetra- 
hedral symmetry (T 5 ) would occur at the softening pres- 
sure of the M§ phonon. This finding is at variance with 
experimental data which indicate that the high-pressure 
phase has an orthorhombic symmetry, compatible with a 
distortion along the PI line. 

The above considerations hold in the hypothesis that 
the strain state of the crystal is constant across the tran- 
sition, except for the isotropic compression due to the 
application of a hydrostatic pressure. We have seen, how- 
ever, that the softening of the Mr" phonon mode is closely 
related to the softening of the shear constant of the crys- 
tal, so that a strong coupling between the soft mode and 
macroscopic strain (i.e., between zone-center and zone- 
border acoustic phonons) is to be expected. Let us con- 
sider now the expression of the Landau crystal energy up 
to fourth order in the phonon amplitude, including the 
coupling with the anisotropic components of the strain 
tensor, e: 

T(u,e) = ^(u,e = 0) + Q(u,e), (3) 
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where Q(u, e) is a polynomial in u and the components 
of the strain tensor containing all the possible products 
up to the fourth order: e 2 , u 2 e, u 2 e 2 , e 3 , and e 4 . All odd 
terms in u must vanish by symmetry. Up to fourth order 
only u 2 e must be considered; in fact the actual value of 
e is the one which minimizes Eq. (|j|) for a given value of 
u, i.e. e oc u 2 so that u 2 e 2 and e 3 are actually of sixth 
order in u while e 4 is of eighth order. To find the form of 
the u 2 e invariants we have first determined all the second 
order monomials, UiUj, such that k^ + kj = G. There are 
nine such monomials which realize a reducible represen- 
tation of the Oh point group. The decomposition of this 
reducible representation into irreducible components is 



displayed in Table VI. The six independent components 
of the strain tensor split into three irreducible represen- 



tation of the cubic group, as indicated in Tabic VII. The 
desired invariants are finally obtained by coupling each ir- 



reducible representation from the uu set (Table VI), with 
an equal irreducible representation from the e set (Table 



VII), and listed in Table VIII. The resulting expression 
for the crystal energy including the coupling with the 
strain reads: 

6 

F(u, e) = + A 2 ^ u i + E Ci o ei<i o 

i— 1,6 i ,j= 1 

+ A ± P 4(u) + E S ^(u, e) + 0(u 6 ), (4) 

i=l,4 i=l,4 

where c,j are the elastic constants, 6j components of the 
strain tensor, and the Q\ are the polynomials listed in 



Table VIII 



Numerical values of the B coefficients have been de- 
termined along lines similar to those used to determine 
the Aj's. The Pll minimum previously found to be sta- 
ble ignoring the phonon-coupling strain is unaffected by 
such a coupling, while the minima along four out of the 
remaining six directions of Table ^ are slightly modified, 
still maintaining the saddle-point character they had in 
absence of such a coupling. The coupling to the strain 
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changes instead the character of the directional minima 
along PI and P2. Let us consider the simpler case of the 
PI minimum. Projecting Eq. (^) along PI one obtains: 

F pl {u, e a , e 4 ) = A 2 u 2 + A P1 u 4 + -c s e 2 s + -c^el + 

(B^es+Bfh^ + Oiu 6 ), (5) 

where u — u\ is the amplitude of the phonon displace- 
ment along (110), B P1 = \B X - 2B 2 , B P1 = B 4 , 
£ s = \{<^xx + £ yy - 2e zz ), and e 4 = e xy . By eliminat- 
ing ei and 62 from Eq. ([5]) by the equilibrium condition 
dJ- pl /de s = dT pl /de^ = 0, one obtains: 

( s = ; £4 = , (6) 

C s C44 



T pl {u) = A 2 u 2 



gPi 2 B pl2S 



Eq. (Q) shows that the coupling between the soft phonon 
and macroscopic strain renormalizes the fourth-order co- 
efficient, making it large and negative, whenever c s or C44 
are small enough. Due to the ongoing softening of the 
shear constant, c s , we find that at the softening pressure 
the fourth-order coefficient is negative, and we conclude 
that the transition must then be first-order, occurring at 
a somewhat lower pressure (see the discussion at the end 
of Sec. HA). A similar behaviour is observed along P2. 
A thorough study of the transition might be performed 
by considering the expansion of the crystal energy up to 
sixth order in the order parameter, and by fitting the 
relevant coefficients to first-principles calculations, as it 
was done when the coupling to the strain was neglected. 
However the presence of many sixth-order invariants and 
the smallness of the computed energy differences make 
such a procedure unpractical. We have preferred instead 
to concentrate on the PI line which is geometrically sim- 
pler and whose symmetry is compatible with the observed 
x-ray diffraction pattern of the high-pressure phase. For 
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this direction we have performed straight energy mini- 
mizations with respect to u and e. By repeating such a 
minimization for different volumes, one directly obtains 
the equation of state of the crystal in the low-symmetry 
phase. 

3. Relative stability 

Csl and CsBr both display a softening of the M§ 
phonon frequency and the ferroelastic instability respon- 
sible for the cubic-to-tetragonal transition, whereas CsCl 
does not show any tendency to phonon softening and 
a very weak tendency towards the ferroelastic instabil- 
ity. We conclude that the high-pressure stable phase of 
CsCl in the pressure range up to « 100 GPa is cubic B2, 
whereas the relative stability of the tetragonal and or- 
thorhombic phases must be checked for CsBr. 

We have performed complete energy minimizations 
with respect to atomic displacements and strain parame- 
ters for both the tetragonal and the newly determined or- 
thorhombic structure. The structural parameters of the 
orthorombic phase of Csl as functions of pressure/ volume 
are shown in Fig. |?]. All these, including volume, are dis- 
continuous at the transition pressure as a consequence 
of its first-order character. The volume discontinuity is 
however very weak (« 0.1%) and it is not visible on the 
scale of the figure. Note also the different magnitude of 
Cesium and Iodine displacements and the weak depen- 
dence oib/a upon applied pressure. The different behav- 
ior of c/a — 1 and b/a — 1 with respect to pressure is due 
to the fact that the former is inversely proportional to 
the nearly vanishing shear constant, c s , whereas the lat- 
ter depends on C44 which is regular in this pressure range. 
We remark that, due to the different amplitude of Cs and 
I displacements in the unit cell, the space group of our 
final structure is different from the G\ v originally 

proposed by Mao et at Q. However, the diffraction pat- 
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tern corresponding to the two groups is the same, and our 
proposed structure is compatible with the experimental 
data. 

We have compared the enthalpies of the various 
phases and identified the transition pressures through 
the Maxwell construction. Results of this procedure are 
shown in Fig. || where the enthalpies of Csl and CsBr 
are shown for the tetragonal and orthorhombic phases 
relative to the cubic one. In agreement with experimen- 
tal findings the high-pressure stable phase of Csl is the 
orthorhombic one (for pressures larger than 21 GPa). 
At variance with the case of Csl, in CsBr the tetragonal 
structure turns out to be more stable all over the explored 
pressure range. The phonon softening is preceeded by 
the cubic-to-tetragonal transformation at Ptet = 58 GPa, 
and we have verified that the tetragonal phase is al- 
ways favored. So we can definitely assign the high- 
pressure phase of CsBr to be tetragonal. The cubic-to- 
tctragonal transition has been experimentally observed 
in CsBr at 53 ± 2 GPa corresponding to a residual vol- 
ume v te t ~ 0.54. H As we have seen, CsCl does not 
manifest any phonon softening upon increasing pressure. 
This compound shows only a weak tendency to a tetrago- 
nal instability at Ptet ~ 69 GPa, but the tetragonal phase 
seems not to be stable at higher pressures. Our results 
are in agreement with this measurement. In CsCl, a weak 
tetragonal instability has been observed experimentally 
at 65 ± 5 GPa, but an extensive study at higher pressure 
is still lacking. 

IV. CONCLUSIONS 

In this paper we have shown that a combined use of 
the Landau theory of phase transitions and accurate first- 
principles calculations based of DFT-LDA provide a re- 
liable scheme for predicting the structural properties of 
ionic systems at high pressures. In the specific case of 
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Cesium halides, we have observed a trend in the high- 
pressure phases which can be related with the polariz- 
ability of the anion. Our results indicate that the more 
polarizable is the anion the larger is the tendency to- 
wards a tetragonal instability of the cubic B2 structure 
(passing from CsCl to CsBr) and, for even larger polariz- 
ability towards an orthorombic structure (from CsBr to 
Csl). Some investigations which we have performed also 
for other alkali halides indicate that when the cation is 
lighter (and hence less polarizable) than Cs, the stable 
high-pressure phase seems always to be the cubic B2 one. 
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TABLE I. Comparison between theoretical and experimen- 
tal lattice constant, ao, bulk modulus, Bo, and electronic 
static dielectric constant, eoo. 



iii, 



Ik] 



Bp [GPa] 



Csl 



theory 
expt. 



AAA 
4.56 a 



15.4 
13.5 a 



3.68 
3.02 c 



CsBr 



theory 
expt. 



A.17 

4.28 b 



20.7 
17.9 b 



3.36 
2.83 c 



CsCl 



theory 
expt. 



3.99 
4.12 c 



24.3 
22. 9 d 



3.16 
2.67 c 



*Ref. 



D Ref. 
c Ref. 
d Ref. 
c Ref. pi 



11 
~| 

ft 
281 



TABLE II. Comparison between calculated and experimentally observed phonon frequencies of Cesium halides at zero 
pressure, at the F, M, X, and R points of the BZ. Units are THz. The uncertainty on the last figure in the experimental data 
is indicated in parenthesis. 





^To(F) LUlo(T) UJTA(M) 


ula(M) 


u>to(M) 


lolo(M) 


ljta(X) 


cjla(X) loto(X) 


LOLO(X) 


ujla(R) 


LOLO(R) 


Csl 


theory 2.09 2.85 1.17 


1.30 


1.48 


2.35 


1.38 


1.59 2.18 


2.21 


1.71 


2.00 




expt. 1.91(5) a 2.74(l) b 1.29(6) a 


1.34(5) a 




2.26(5) a 


1.21(5) a 


1.30(6) a 2.07(5) a 


2.16(6) a 


1.75(5) a 




CsBr 


theory 2.44 3.53 1.27 


1.48 


1.98 


2.79 


1.47 


2.15 2.42 


3.00 


1.89 


2.63 




expt. 2.29(2) c 1.22(3) c 


1.61(3) c 


1.77(5) c 


2.69(5) c 


1.32 c 


1.91° 




1.89 c 


2.40 c 


CsCl 


theory 3.33 5.06 1.24 


1.86 


3.18 


4.04 


1.43 


2.44 3.21 


4.50 


1.91 


3.94 




expt. 3.17(2) d 1.19(3) d 


1.90(2) 


2.77(2) d 


3.75(3) d 


1.31 d 


2.98 d 




2.02 d 


3.85 d 



a Ref. 
b Ref. 



30 
321 



c Ref. [!3| 
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TABLE III. Phonon wavevectors, q, and polarizations, e, 
of the six components of the order parameter, u. 

Ml U2 Us U4 Us, «6 

q (oil) (oil) (p|5 (pi) (ip) (ipT 

e ( 10 °) (°7f^) (° 10 ) (ti ^) (° 01 ) (t5^ 



TABLE IV. Fourth-order invariant polynomials which en- 
ter the free energy expansion, Eq. ^. 

Pl(u) uf + u\ + uj + ui + uf, + ut 

r>2 / \ 22,22,22 
P 4 (u) U±U 2 + U 3 U 4 + U 5 U 6 

P|( u ) — uiu 2 (ul + u\ — u\ — u|) + 

uzUa{u\ + u\ — u'{ — ul) + 

/ 2 i 2 2 2\ 
tt5«6(Wl + U 2 — «3 — U 4 ) 

Pi (u) — U1U2U3U4 — M1M2W5M6 + U3U4U5U6 
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TABLE V. Classification of the possible stable phases 
along particular directions of the order-parameter space. 



Direction Space group 

PI ( 1, 0,0, 0, 0, 0) Dfh 

P2 ( 1,-1,0, 0, 0, 0) D\l 

P6 ( 1, 0,1, 0, 1, 0) Cl 

P7 ( 0,-1,0, 1, 0, 1) D 7 3 

P9 ( 1,-1,1,-1, 0, 0) Dll 

P10 ( 1,-1,0, 0, 1,-1) Dll 

Pll ( 1,-1,1, 1, 1, 1) T 5 



TABLE VI. Decomposition of the set of lattice-periodic 
quadratic monomials in the order parameter into irreducible 
representations of the cubic group. 



Label Basis functions 



—U\U2 + U3U4 + U 5 U(j 

T~3 — 2u\u z — U3U4 — u^ue; 

W3W4 — u$ue 

r+ u\ + U% + U3 + u\ + u\ + u\ 

r+ 2(u\ + ul) - ul - u\ - ul - u%- 

2,2 2 2 

u 3 + w 4 - u 5 - u 6 

r+ 2 2 2 2 2 2 

J u- L - u 2 ; u 3 - u 4 ; u 5 - u 6 



TABLE VII. Decomposition of the set of independent com- 
ponents of a symmetric rank-2 tensor into irreducible repre- 
sentations of the cubic group. 

Label Basis functions 

£xx ^yy ) 2a Z z £xx £m) 

£xy] €yz\ €zx 



TABLE VIII. Fourth-order polynomial invariants of the 
phonon-strain coupling. 

Ql(u,e) (e X x + £yy + £zz){u\ + u{ + ui + u\ + U% + U%) 

Qt{u,e) (2(u\ +ul)-ul-ul-ul-u%)x 

(2e Z z — (-xx ~ tyy) + 3(ll 3 + U 4 — U 5 — U 6 )(e xx — £yy) 

Ql(u, e) ( — 2uiU2 — u 3 U4 — ur,u G )(exx — £ yy ) + 

(-'U3U4 + u 5 « 6 )(2e zz — e xx - <Lyy) 

Qt{u, e) (u\ - u\)t xy + [u% - ul)e yz + («§ - u%)e xz 
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FIG. 1. Sketch of the Csl cell with the atomic displace- 
ments corresponding to an acoustic phonon. The magni- 
tude of Cesium displacements with respect to Iodine displace- 
ments is exaggerated for clarity. The solid lines indicate the 
cell of the distorted (orthorhombic) structure, while the cell 
of the undistorted (cubic) structure is indicated by a dashed 
line. The shaded area correspond to the (110) plane. 

FIG. 2. Phonon dispersion relations along some symmetry 
directions of the Brillouin zone in Csl at equilibriumi volume. 
Experimental data indicated by open and black circles are 
from Ref. j3o| ; an infrared measurement for lulo(T) is indi- 
cated by an open square. [^2) 

FIG. 3. Phonon dispersion relations along the E (110) 
for Csl at equilibrium volume and just above and below the 
softening pressure of the acoustic mode. 'Negative' fre- 
quencies actually mean 'imaginary' (i.e. negative squared fre- 
quencies) . 



FIG. 4. Frequencies of the acoustic phonon as func- 
tions of the molar volume for Csl, CsBr and CsCl. Arrows 
indicate the softening volume. P* is the softening pressure. 

FIG. 5. Enthalpies vs c/a for the three compounds corre- 
sponding to a pressure Pi < P te f, P% ~ Ptet; Pz > Ptet- 

FIG. 6. c/a parameter as a function of pressure for Csl and 
CsBr. Ptet indicates the transition pressure to the tetragonal 
phase. 

FIG. 7. Structural parameters of the orthorhombic phase 
of Csl as functions of different values of the volume/pressure. 
Arrows indicate the transition pressure/volume. The vol- 
ume discontinuity at the transition is not visible on the scale 
of the picture. P ort indicates the transition pressure, while 
v cu b and v or t are respectively the volumes of the hig h- and 
low-symmetry phases at the transition (from Ref. |12|). 




FIG. 8. Enthalpies of the orthorhombic and tetragonal 
phases, relative to the enthalpy of the cubic phase, as func- 
tions of the applied pressure in Csl and CsBr. P ort and Ptet 
are the transition pressures from the cubic to the tetragonal 
and to the orthorhombic phases respectively. 
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